Universality of the metal-insulator transition in three-dimensional disordered systems 



O 
O 

(N 

(N 



C 
i 

C/3 



I 

c 

o 
o 



(N 
> 

in 
o 

O 
O 

03 



o 
o 



X 
S3 



J. Brndiar and P. Markos 
Institute of Physics, Slovak Academy of Sciences, 845 11 Bratislava, Slovakia 

The universality of the metal-insulator transition in three-dimensional disordered system is con- 
firmed by numerical analysis of the scaling properties of the electronic wave functions. We prove 
that the critical exponent v and the multifractal dimensions d q are independent on the microscopic 
definition of the disorder and universal along the critical line which separates the metallic and the 
insulating regime. 

PACS numbers: 73.23.-b, 71.30., 72.10. -d 



One of the main problem of the disorder induced metal- 
insulator transition (MIT) is the proof of its universal- 
ity. In the pioneering work 0, it was conjectured that 
if the sample size exceeds all the length parameters of 
the model, then the conductance, g, is the only parame- 
ter which governs MIT. This scaling hypothesis has been 
confirmed by various numerical analysis, with the help of 
the finite-size scaling 0, 0] ■ 

Generally accepted scenario of the Anderson localiza- 
tion is that disorder broadens the conductance band. 
Electron states in the tail of the band become localized, 
separated from delocalized (metallic) states by the mo- 
bility edges, E c . System exhibits the MIT if the Fermi 
energy, Ep, crosses the mobility edge. With increased 
disorder, E c moves towards the band center. There is a 
critical value of the disorder, W c , for which E c reaches 
band center, E c = 0. For disorder W > W c , all electronic 
states inside the band become localized. Phase diagram 
in the energy-disorder plane was calculated in |4j and is 
schematically shown in the upper panel of Fig. ^ 

At the band center, E = 0, the universality of the MIT 
was confirmed by detailed numerical analysis of the disor- 
der and system size dependence of Lyapunov exponents 
in quasi-one dimensional systems J3l |5j, mean conduc- 
tance 0] , conductance distribution jj , and level statistics 
0, . These studies determined the value of the critical 
exponent v, which determines the divergence of the cor- 
relation length, i ~ \W-W C \-", asv = 1.57±0.02 5, 6]. 
The analysis of MIT along the critical line (non-zero en- 
ergy E) is more difficult because the critical region is 
narrower and finite size effects are stronger 01 • Critical 
exponent, v, was obtained only in models with random 
hopping , and very recently in [l^ . 

In this paper, we present numerical proof of the univer- 
sality of MIT. By scaling analysis of the electronic wave 
functions in the vicinity of two critical points, shown in 
the upper panel of Fig. ^ we prove that the critical ex- 
ponent v and fractal dimensions d q of the wave function 
(defined below) are universal along the critical line. 

Electron eigenenergies and wave functions are calcu- 
lated for three-dimensional Anderson Hamiltonian. 



Here, r counts the sites of the three-dimensional (3D) 
lattice of the size L 3 , e r is the random energy dis- 
tributed either with the Box distribution, -Ps(e) = 
(2/W)9(W/2 - |e|) or with the Gaussian distribution, 
P G (e) = V2ttW 2 exp ~(e 2 /2W 2 ). Parameter W mea- 
sures the strength of the disorder. For E = 0, the critical 
disorder W c « 16.5 (6.15) for the Box (Gauss) distribu- 
tion of random energies, respectively. 

The quantities of interest are inverse participation ra- 
tios (IPR), Ig(E n ). By definition, [3 
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FIG. 1: (Color online) :w Top panel shows schematic phase 
diagram for 3D Anderson model. Solid line separates metal- 
lic (shaded area) and localized states. Open circles shows the 
position of critical points studied in this paper. Note that the 
mobility edge, E c lies outside the unperturbed energy band. 
Two bottom panels present numerical data for I2, given by 
Eq. for cubes of the size L — 8 (middle) and L = 16 (bot- 
tom) and with Gaussian disorder W — 2 = 0.325 W c . In the 
metallic region, I2 decreases when L increases, while I2 be- 
comes i-independent in the tail of the band, where electronic 
states are localized. 
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Here, E n and $ n ( r ) is the nth. eigenenergy and eigen- 
function of the Hamiltonian respectively. 

The size dependence of IPR indicates the character 
of the eigenstate. If the nth eigenstate is conductive, 
the wave function is distributed throughout the sample 
and |$ n (r)| oc L~ d / 2 . Inserting in Eq. J3J) we obtain 
that I q (E n ) oc L^ 1 -^. For localized state, $ n (r) is non- 
zero only in a small region, where |$„(r)| ~ 1. Hence, 
I q (E n ) ~ 1, too. These size dependences are shown in 
two bottom panels of Fig. ^ 

Different size dependence of I q (E) enables us to use 
IPR as a scaling variable for the calculation of the critical 
parameters in the same way as Lyapunov exponents Q, 
the conductance Q, or level statistics ||. However, in 
contrast to the Lyapunov exponent or mean conductance, 
IPR is not a size inde pen dent constant at the critical 
point, but decreases as [Mllallq 
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where d q are fractal dimensions. This makes the scal- 
ing analysis slightly more difficult. On the other hand, 
fractal dimensions, d q , represent a new set of parameters, 
which can be used for the verification of the universality 
of the MIT. We expect that d q are universal constants 
for all critical points along the critical line. 

The energy spectrum of the Hamiltonian depends on 
the system size, L, and on the microscopic details of the 
disorder in a given sample. Therefore, we have to calcu- 
late an average values, defined as follows. For each sys- 
tem size, we consider a statistical ensemble of N s samples 
which differ only in the realization of random energies, e r . 
For each sample, we calculate all eigenenergies, E n , lying 
in a narrow energy interval, E — S, E + S, and calculate 
corresponding I q {E n ). For the ith sample, the number 
of eigenstates, rii, depends on the microscopic realization 
of the disorder. Also, since the values of I q might fluc- 
tuate in many orders in magnitude in the critical region 
[lfij (these fluctuations are shown in Fig. EJ, it is more 
convenient to study logarithm of I q . Thus, our scaling 
variable is then defined as 
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where A^ ta t = J2i n i- in our calculations, AT sta t ~ 10 5 
(5 x 10 4 for L > 50) and S = 0.025. With these param- 
eters, we calculate Y q (E) with relative accuracy better 
than 0.2%. Numerical data were collected by LAPACK 
subroutines for L < 16. For larger system size (L < 54) 
we use our own program based on the Lanczos algorithm. 

We expect that Y q is a good scaling variable, so that 
it behaves in the vicinity of the critical point as 



Y q {E, L) = Y q c - d q InL + A(E - E C )L X I V ', 
for the fixed disorder W, and 
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FIG. 2: (Color online) Ys(E,L) as a function of system size 
L for various energies E — 6.50, 6.55, . . . , 6.80 (from bottom 
to top). Thin solid lines are fits, Eq. (y. Thick solid line 
is (lnl 5 ) = Yg — 4 d 5 lnL (Eq. [7J for the critical energy 
Ec w 6.58. The same fit for other three critical points are 
also shown (described by legend). Data confirm universal 
L-dependence of F5 for three critical points. The different 
behavior for the critical point BE C is due to strong finite-size 
effect, discussed in the text. 



for the fixed energy E = 0. 

For a given q, we fit obtained numerical data for Y q to 
Eqs. © or JBJ. Obtained results are in Figs. [3 El and E| 
Typical data for q = 2 and q = 5 are given in Table [I] 

Figure shows the L-dependence of Y5 for Gaussian 
disorder W = 2 and for energies close to the critical 
energy, E c 1=3 6.58. We see that at the critical point, 
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TABLE I: Critical exponent, v, and fractal dimensions, d q , 
calculated by scaling analysis of the inverse participation ra- 
tio, I q , for four critical points. Calculated position of critical 
points, E c , W c , is given in the 3rd column. GW C - Gaussian 
disorder, band center, GE C - Gaussian disorder, band tail, 
BWc - Box disorder, band center, BE C - Box disorder, band 
tail. Data for system of the size L m i n < L < L max were used 
in the scaling analysis. 
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E = E c , Y q (E c ) decreases logarithmically, 
Y q (E c ,L) = Y q c -d q ]nL, 



(7) 



in agreement with Eq. (JSJ. Outside the critical point, 
the L dependence of Y q changes due to the presence of 
the term ~ (E — E^L 1 ^ . Similar analysis, performed 
for other critical points confirms the universality of the 
relationship (0). This indicates that parameters Y q) d q 
and A are universal. 

There are two sources of inaccuracy of the scaling anal- 
ysis: (i) if the critical region is not sufficiently narrow, 
then the linear term ~ E — E c could not be sufficient to 
describe a correct independence of numerical data and 
higher order terms of the expansion (0 must be consid- 
ered 0| • To test the accuracy of the linear approximation, 
we add also cubic term, ~ (E-E c ) 3 in Eq. ©. We found 
that such correction does not influence obtained critical 
parameters and might be neglected, (ii) for small sys- 
tem size, L, the variable Y q suffers from finite-size effects 
(FSE) jg. The role of FSE can be estimated by the scal- 
ing analysis for data of restricted system size, L > L alln . 
For three critical points, GE C , GW C and BW C (the posi- 
tion of critical points is given in Table QJ, we found that 
data for L > L m ; n ~ 16 — 20 are already free of FSE. 
However, we were not able to obtain reliable critical pa- 
rameters for the critical point BE C . As shown in Table |U 
and Fig. our numerical data for this critical point are 
still far from limiting values. For this critical point we 
can only demonstrate the convergence of critical param- 
eters to expected universal values when L mm increases 
(inset of Fig. ■ Larger samples are necessary to prove 
this convergence numerically 

Figure summarizes our data for fractal dimensions, 
d q . Data confirm that the spatial distribution of the wave 
function is universal, independent on the position of the 
critical point along the critical line. Presented data for 
d q are in agreement with previously reported values . 

Figure 0] presents obtained data for the critical expo- 
nent, v. Wee see that v converges to the generally ac- 
cepted value, v ~ 1.57 when either L m - m or q in- 
creases. In order to show the effect of the system size, we 
plot for each value of q a few data obtained with increas- 
ing minimal system size used in the scaling procedure. 
Inset of Fig. 01 presents our data for the critical point 
BE C . Because of strong finite-size effects, much larger 
systems are necessary for the estimation of reliable val- 
ues the critical exponent. 

To understand the origin of the FSE, we calculated 
the density of states, p(E), in the critical region of all 
four critical points. As shown in left Fig. [S] p(E = 0) 
changes only in a few % when disorder varies around 
the critical value, W c . Contrary to the band center, the 
density p(E) around the mobility edge GE C decreases sig- 
nificantly (almost by factor of two) in the critical region 
(right Fig. |3J. In the case of critical point BE C , this de- 
crease is even more significant. Intuitively, one expects 
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FIG. 3: (Color online) Fractal dimensions, d q , calculated by 
scaling analysis of Y q for three critical points. A small hori- 
zontal shift of symbols for a given q indicates minimal system 
size, L m in, used in the scaling analysis: Z/ m i n = 8 (left), 16 
(middle) and 20 (right). 
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FIG. 4: (Color online) Critical exponent, u, calculated from 
the scaling behavior of Y q for three critical points. A small 
horizontal shift of symbols for a given q indicates minimal 
system size, L m j n , used in the scaling analysis: L m i n = 8 
(left), 16 (middle) and 20 (right). Inset shows the critical 
exponent v calculated for critical point BE C from Y2 (circles) 
and Ys (triangles), v increases when data for smaller system 
size are omitted (L m m increases). This indicates that the 
deviation from universality, observed in this critical point, is 
only finite size effect. Solid (dashed) lines are fits ao + ai/L a2 
with ao Ri 1.36 (1.29) for q = 5 (g = 2), respectively. 



that the one-parameter scaling works better when elec- 
tronic states inside the critical region have the same, or 
at least comparable, density. Since the interval of ener- 
gies, AE = 0.3, used in the case of GE C is already suf- 
ficient to get correct critical parameters, we expect that 
the the energy interval for BE C must be AE w 0.2. How- 
ever, narrower interval of energies requires larger system 
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FIG. 5: The density of states, p(E), arround the critical 
points. Dashed vertical lines indicates the position of critical 
points. Left: the disorder dependence of the density p(E = 0) 
at the band center. p(E = 0) changes only in 10% when dis- 
order varies wihin the intervals (15.5 < W < 17.5 for BW C 
and 6.5 < W < 6.8 for GW C ) used our scaling analysis. Left: 
The density of states, p(E) around the critical points GE C 
and WE C . Solid lines are linear fits with slopes 0.018 (GE C ) 
and 0.026 (BE C ). 
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FIG. 6: Probability distribution P(ln/5) for three critical 
points and L — 40. We used data from narrow energy in- 
terval, E c ± 5E. The size of the system is L = 40. Note 
that neither mean value nor the distribution is system-size 
independent. Mean value, (ln/5) = —13.08. The value of the 
variance, var In J5 = 11.6, agrees with estimation of Ref. [lr^|. 
Solid line is Gaussian distribution with the same mean value 
and variance. 



size. To have the term (E — E C ]L 1 I V of the same magni- 
tude as in the GE c -ca.se, the size of the system must be 
(3/2)i/ ~ 1.8-times larger than in the GE C case. Conse- 
quently, we expect that the scaling analysis of systems of 
the size L max ~ 90—100 would lead to correct estimation 
of the critical exponent also for the critical point BE C . 

We expect that I q is less sensitive to finite size effects 
when q increases. Since the spatial distribution of elec- 
tron is non-homogeneous in the critical regime, the main 
contribution to I q is given by sites with large values of 
| $„ | . Those sites are well localized in space and so insen- 
sitive to the system size. This is confirmed also by data 
shown in Fig. 0] 



The universality of the spatial distribution of the crit- 
ical wave function is confirmed also by Fig. [B] which 
shows the probability distribution of In I5 for three crit- 
ical points, GE C , GW C and BW C . The width of the dis- 
tribution, a q = var In I a ~ 3.4, is close to the limiting 
value reported in Ref. 17] . However, the distribution 
P(ln/ g ) is not system size invariant. The mean value 
of ln/q decreases when system size increases, while the 
distribution always possesses a long tail for lnl q ~ 1 

In conclusion, we investigated numerically the wave 
function of electron in the critical regime of the metal- 
insulator transition. By the scaling analysis of the log- 
arithm of the inverse participation ratio, I q , we calcu- 
lated the critical exponent, v, and fractal dimensions, 
d q . We found that these parameters depend neither on 
the microscopic details of the disorder nor on the posi- 
tion of the critical point. This result confirms that the 
metal-insulator transition is universal along the critical 
line which separates metallic and insulating regimes. 

This work was supported by grant APVV, project n. 
51-003505 and VEGA, project n. 2/6069/26. 
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